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Abstract 


It is shown that the effective 4D lattice Dirac operator of optimal 
lattice domain-wall fermions with finite Ng (in the fifth dimension) is 
exponentially local for sufficiently smooth gauge background. 
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Recently, it has been shown that one can construct a lattice domain-wall 
fermion action such that the effective 4D lattice Dirac operator preserves the 
chiral symmetry optimally for any given hnite Ng in the hfth dimension |^. 
Explicitly, the action of optimal lattice domain-wall fermions readsQ 
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= + {uJsD^{x,A) - l)P_6s',s+l 

s,s'=l x,x' 
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+ {u;sDig{x,x') - 1) P+6 s’,s-ili’ix', s') 
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( 1 ) 


with boundary conditions 

771 771 

P+ij{x,0) = -^P^ij{x,Ns), P-'iP{x,Ns + 1) = -^P.^{xA) , 

ZtTIq ZTTT/Q 

and the quark fields with bare mass rUg can be constructed from the left and 
right boundary modes 


q{x) = 1) -h P+'ipix, Ns) 

q{x) = 'tl){xA)P++ i’{x-,Ns)P- . 

Here is the 4D Wilson-Dirac operator with a negative parameter (—mo), 

4 

Dw = + W -mo, mo e (0, 2) 

t^(,X,X ) —)U^(^x)6x',x+fi Pjii.X )6x',x—fi\ 

W(x, x') = ^ i [26x,x' - U^{x)6x',x+t, - Ul{x')6x',x-^J] , 

M=i ^ 

and the weights 

x’s = T—- P^sn^ (n^; N), s=l,---,Ns (2) 

'^min 

where sn(ns; k') is the Jacobian elliptic function with argument Vg and mod¬ 
ulus P = — AP„/ \max) i^max ^min ^lie maximum and the minimum 

of the eigenvalues of H^, Hyj = The weights {a;^} are obtained from 

the roots {us = s = 1, ■ • •, Ng) of the equation 

J ^ “ ^R^z’'^\u) = 0 , Ng = 2n + l 

^ \ I - ^R^z~^''''\u) = ^ , Ng = 2n 

Hn this paper, we suppress the lattice spacings (a and 05), as well as the Dirac and color 
indices, which can be restored easily. 
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The argument Vg in is 


V 


S 


where 


sn 


-1 


' 1 + 3A 

(ITaP 


;VT^] + 


2K' 



(3) 

(4) 

(5) 


K' is the complete elliptic integral of the first kind with modulus k\ and 0 is 
the elliptic theta function. Note that since sn^(;) < 1. 

After adding pseudofermions (Pauli-Villars helds) with hxed mass = 
2mo to the action one can derive the effective 4D lattice Dirac operator 
for the internal quark loops 


where 


D{mg) =r{D^ + mg){l + rDyj \ 


1 

r =- 

2mo 


_ 1 'y^^opt 

1 'ys^opt 

I H^Ry^\Hl), Ng = 2n + l 

\ Ng = 2n, 


( 6 ) 


(7) 


and Rz{H^) is the Zolotarev optimal rational approximation for the inverse 
square root of H^,, 


Tj(n,n) 

Hz 


{Hi 


dp -pr 1 + hi,/C 2 I 

^min [—1 1 + hllc2i-i 


( 8 ) 


and 


p{n-l,n) 

Hz 


{HD 


d', Yippi + hi/d,;) 
x^in nr=i(i + ^y4-i) 


(9) 


in which = Hl^/XH^, and the coefficients do, q and c'l are expressed in 
terms of elliptic functions 0 depending on Ng and XDx/Xmin- 

In the limit Ng —> 00 , Sopt —> H^(HD~^^^, thus rDc (0) becomes chirally 
symmetric. Therefore, in the massless limit {irtg = 0) and Ng —>■ cx), D is 
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exactly equal to the overlap Dirac operator |^, and satisfies the Ginsparg- 
Wilson relation p 


-D75 + 75 -D = 2D75D . 


This implies that the effective 4D Dirac operator of optimal DWF is expo¬ 
nentially local for sufficiently smooth gauge background, and is topologically- 
proper (i.e., with the correct index and axial anomaly), similar to the case of 
overlap Dirac operator. At finite Ng, the effective 4D lattice Dirac operator 
is exactly equal to the overlap Dirac operator with approximated by 

Zolotarev optimal rational polynomial. 

The valence quark propagator coupling to physical hadrons is 

{Dc + mq)~^ = r(l - rmq)~^[D~^{mq) - 1] , 

where D~^{mq) can be computed via the 5-dimensional lattice Dirac operator 
of optimal DWF. 

In this paper, we show that the effective 4D lattice Dirac operator (|^) is 
exponentially local for sufficiently smooth gauge configurations. Since Eq. (H) 
can be rewritten as 


D{mq) = rniq + (1 — rmq)D , 


where 

D = rDc(l -h rDc)~^ = ^(1 -h i^Sopt) , 

thus it suffices to assert the exponential locality of D ([TlD . 
First, we rewrite (H) and @ in partial fractions 


where 




S TT 75777 

^ bi 

hi 1 + hl/c2i-i 


Cl 


hi 


l - sn 2(^;«0 ’ 

J_2A 1 nr=~/(l-C2;-i/c20 

Amm 1 -|- A M nr=l,i^z(l ~ C2l-l/c2i-l) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


A and M are defined in (^) and (^ respectively. 

Therefore, it suffices to examine the locality of each term in (^) and (^), 
by expanding each fraction in terms of orthogonal polynomials, similar to 
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the approach used in 0. Here I expand each fraction in term of Chebycheff 
polynomials of the hrst kind 


2 C 2/-1 


l + ( 6 -l)sinh 6 'z 


1 + 2 ^ 

^ k=l 


6 -"^' T,{z) 


where 


\2 

^ _ ^max 




d] = cosh 


-1 ( b + 1 + 2c2z_i 


6-1 


2 ; = 


6 + 1 - 

6-1 


(16) 

(17) 

(18) 
(19) 


The expansion (|l 6 l) is valid provided that 6 *; 7 ^ 0 (i.e., > 0), which is 

ensured by any gauge conhguration satisfying 0 


plaquette ||l-t/(p)||<6< 


1 — |1 — mol 

6(2 + 72 ) 


( 20 ) 


since 

(v^l - 6(2 + 72)e - |1 - m„|D" < AV„ < < (|4 - m„| + 4)" . (21) 


Because only has nearest neighbor coupling, it follows that h^(a;, ?/) = 0 
for all separations bigger than two lattice spacings, i.e., \x^ — yy\ > 2a. 

Thus, for the Chebycheff polynomial Tk{z) of degree k, [Tk{z)]{x,y) = 0 for 
T,p \Xfi - Vfil > ‘^ka. Then, (|T|) gives 


1 + 6 ,^/ C 21-1 


{x,y) = 


2 C 2 i-l 


(6 - 1 ) sinh Oi 


[T,{z)]{x,y) , (22) 


where m = Yp \xp~yp\/‘^o, > 1. Using the triangle inequality and the property 
||7fc(^)|| < 1, we turn Eq. (p2D into the inequality 


1 + hl/c2i-i 


{x,y) 


< 


2 C 2 i-l 


n-me. 


(6 — 1) sinh6*i (1 — e~^‘) 


and it follows that 


< 


< 


2 biC2i-i 1 f sx i 

—1 g smh«,(i-«-“') ? 1^" 

2 biC2i-i 1 \ / 6*1 

6—1 sinh 9i (1 — e~^‘) J \ 2a ^ 



(23) 
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where the last inequality is due to the property { 6 i < 62 <■■■ < On) from 
([ISD, since (0 < ci < C3 < ■ • • < C2n-i) from From (|TT| ) - ([T^) , and it 
follows that D (O) is exponentially-local, and so is D{mq) (pUf). 

Obviously, the above analysis also ensures that the overlap Dirac opera¬ 
tor with approximated by Zolotarev optimal rational polynomial is 

exponentially-local. Note that, in practice, the least upper bound of (|^) is 
much smaller than that on the r.h.s. of (^3]). This implies that D{x,y) is still 
exponentially local even for rough gauge conhgurations bounded by an e much 
larger than its least upper bound in (|^) . 

Further, it is straightforward to assert the differentibility of D o with 
respect to the gauge held, as well as the exponential locality of D in the 
presence of an isolated near-zero modes of H^, thus the details can be omitted. 

Finally, it is instructive to observe that the decay constant 61 in (p3D is a 
monotonic decreasing function of Ng. In the limit Ng —> 00 , Ci —>■ 0, and 61 
becomes the minimum. 


9 = cosh' 


'b + V 

6-1 


= cosh 



In other words, at any hnite Ng, D{x, y) is more localized than itself at Ng = (yo. 
In the limit a;* = 1, the effective 4D lattice Dirac operator is equivalent to 
the overlap Dirac operator with polar approximation, and ci is replaced with 
di = tan‘^{ 7 i/ 2 Ns), thus its decay constant ( 6 * 1 ) is always smaller (i.e., less 
localized) than that of the optimal DWF, for any hnite Ng. 

For the conventional DWF with open boundary conditions, it has been 
estimated |p that the ehective 4D lattice Dirac operator at = cxo behaves 
as 


DeQ{x,y)\ \ < 





However, for hnite Ng, the decay constant of Dgg(a;, y) has not been obtained. 
Note that the factor \x^ ~ Ufj-l) in D^Q{x,y) makes it decay slower than 
purely exponential decay. Nevertheless, this is only an artifact due to the 
expansion [§ in terms of Chebycheh polynomials of the second kind. If one 
expands the fraction in terms of Chebycheh polynomials of the hrst kind as 
(pTl) , then one obtains a purely exponential bound. 


DQs{x,y)\ \ < 
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